Geometric approach to the dynamic glass transition 
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We numerically study the potential energy landscape of a fragile glassy system and find that 
the dynamic crossover corresponding to the glass transition is actually the effect of an underlying 
geometric transition caused by the vanishing of the instability index of saddle points of the potential 
energy. Furthermore, we show that the potential energy barriers connecting local glassy minima 
increase with decreasing energy of the minima, and we relate this behaviour to the fragility of 
the system. Finally, we analyze the real space structure of activated processes by studying the 
distribution of particle displacements for local minima connected by simple saddles. 



Despite a large number of investigations, there is still 
much to understand about the dynamic glass transition 
in supercooled liquids. The basic problem is that, strictly 
speaking, there is no dynamic transition at all. In sys- 
tems known as fragile liquids Q , experiment finds a sharp 
rise of the viscosity in a very narrow interval of temper- 
ature upon cooling. The shear relaxation time increases 
by several orders of magnitude in a few degrees, and it be- 
comes impossible to perform an equilibrium experiment. 
Nevertheless, sharp as this behaviour may be, it is not a 
genuine dynamic singularity. At the other extreme of the 
experimental spectrum we find strong liquids m, which 
experience a gentle increase of the relaxation time, often 
according to the Arrhenius law. Even in such systems 
though, when the viscosity becomes too large, equilib- 
rium can no longer be achieved in experimental times. 

The glass transition temperature Tg is conventionally 
defined as that where the value of the viscosity is lO^'^ 
poise. Below Tg equilibrium experiments become really 
hard to perform and a sample can be considered to be in 
its glass phase. However, Tg is just a conventional exper- 
imental temperature, defined out of the need to mark the 
onset of glassy dynamics. The attempt to give a theo- 
retical description of such an ill-defined "transition" may 
therefore seem pointless. 

On the one hand, this conclusion is correct for the 
strongest liquids: here nothing peculiar happens close 
to Tg, and the glass transition fully displays its purely 
conventional nature. On the other hand, the most fragile 
systems resist such an objection, simply by virtue of the 
extremely steep increase of relaxation time in a small in- 
terval of temperature around Tg. This fact suggests that 
some kind of new physical mechanism is indeed responsi- 
ble for the onset of the glassy phase in fragile supercooled 
liquids. We share this view, and the aim of this Letter is 
to shed some light on the nature of this mechanism. 

The key idea is that the sharp dynamic crossover ob- 
served in fragile liquids is a consequence of an underly- 
ing topological transition, controlled by energy, rather 
than temperature. More precisely, the existence of an 
energy level where the instability index of the stationary 
points of the potential energy vanishes is responsible for 



a change in the dominant mechanism of diffusion. If en- 
ergy barriers are large, this change in the mechanism of 
diffusion causes the fast increase of the relaxation time. 

This study is part of a more general program aimed 
at explaining glassy dynamics in terms of properties of 
the potential energy landscape. Our method generalizes 
the ideas of Goldstein |g|, and Stilhnger- Weber ||], by 
extending to unstable stationary points the analysis for- 
merly restricted to minima of the potential energy. The 
first steps in this direction have been done in El, build- 
ing on the ideas of [glJg], and more concrete results have 
been recently obtained in MM ■ Further inspiration came 
from the approach of Keyes and coworkers [p[, which 
related diffusion to the stability properties of instanta- 
neous configurations. In the present work we firmly es- 
tablish the connection between topological properties of 
the landscape and fragile glassy dynamics. Furthermore, 
we study the role of potential energy barriers and we an- 
alyze the real space structure of activated processes. 

We consider a soft-sphere binary mixture ||l^ , a fragile 
model glass-former. In addition to capturing the essen- 
tial features of fragile glasses |10|-^^ , this model can be 
thermalized below Tg with the efficient MC algorithm of 
p2[ . Furthermore, previous investigations of the saddle 
points have focused on Lennard- Jones systems, so it is 
useful to look at a broader class of models. Most of our 
data are obtained for iV = 70 particles, but we tested 
our key results for N = 140 as well. We impose periodic 
boundary conditions in d = 3 dimensions. Particles are 
of unit mass and belong to one of two species a = 1,2, 
present in equal amounts and interacting via a potential 
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The radii CTq are fixed by 02! o\ — \.1 and setting the 
effective diameter to unity, that is (2cri)'^ -I- 2(cri -I- 172)'^ + 
(20-2)'^ — 4Zq, where Zq is the unit of length. The density 
is p = N jV = 1 in units of l^^, and we set Boltzmann's 
constant ks — 1- We obtain equilibrium configurations 
at several temperatures by the swap Monte Carlo algo- 
rithm of ll^ . A long-range cut-off at re — VS is imposed. 



However, to find the stationary points we need a poten- 
tial with a continuous second derivative. Thus, instead 
of simply shifting the pair potential by a constant Cy (so 
that Vij{r > Tc) = 0), we use a smooth cut-off, setting 
Vij{r) = Bij{a — r)^ for re < r < a and Vij{r) — for 
r > a, fixing a, Bij and Cij by imposing continuity. 
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FIG. 1. Top: Average instability index density k vs. po- 
tential energy density u of the stationary points. Bottom: 
Temperature vs. equilibrium bare potential energy density, 
Ub = Ucq — 3/2 fcsT. Inset: k(u) on the whole range sampled. 
Open symbols: N = 140, filled symbols: A^ = 70. 



We sample the stationary points of the potential en- 
ergy by quenching the equilibrium MC configurations 
onto saddle points. This is done by numerically solv- 
ing the 3A^ nonlinear equations dV /dri = by means 
of a backtracking Newton method with finite-difference 
approximation to the Jacobian p4| . Once a saddle point 
is found, we measure its potential energy U and its insta- 
bility index K, that is the number of negative eigenvalues 
of the Hessian matrix at the saddle. 

In Fig. |l] (top) we plot the average index density 
k — K/3N as a function of the potential energy den- 
sity u — U/N. As in ^, we find a well defined func- 
tion k{u) which vanishes at a threshold value of the en- 
ergy, Mth- For A'^ = 70, a comparison between the lin- 
ear fit of the data and the last point of the curve gives 
Mth — 1-75 ± 0.01. For A^ = 140, despite worse statistics, 
we find Mth = 1.73 ± 0.01, consistently with the result for 
N = 70. Note that the threshold energy is not the ground 
state of the system, and in fact we found minima down 
to Mo = 1.68 (A^ = 70). The threshold energy marks 
the border between unstable saddle points, dominating 
the landscape above uth, and stable minima, dominant 



for uq < u < uth- Therefore, a topological transition 
takes place at Uth, where the stability properties of the 
landscape change. 

A system confined to the minima-dominated region of 
the landscape, u < Uth, raust resort to barrier hopping 
to diffuse in phase space. It is therefore essential to find 
the temperature Tth below which this confinement takes 
place. To this end we must realize that a system trapped 
in a single potential well has a potential energy density 
equal to the bare energy of the bottom of the well, plus 
a vibrational contribution proportional to /cbT. There- 
fore, it is the bare potential energy of the system which 
we must compare with the threshold energy [pj. We 
can write the bare energy as Ub{T) =^ Ucq{T) — 3/2 ksT. 
When the bare energy Ub{T) drops below the threshold 
Mth, the system is effectively confined to the minima dom- 
inated region of the landscape. For A^ = 70 this happens 



at Tth = 0.242±0.012 (Fig 
previous investigations 



bottom). Note that, unlike 
for A^ = 70 we thermalize 



the system below the threshold, giving an accurate de- 
termination of Tth Jljl • Hence, the dynamic effect of the 
topological transition at Mth must be a qualitative change 
in the mechanism of diffusion at Tth- 

The fact that barrier crossing becomes the main mech- 
anism of diffusion below Tth does not necessarily imply 
a slowing down of the dynamics: large energy barriers 
at Uth are also needed, i.e. AU{uth) substantially larger 
than fcsTth- The value of AC/(uth) can be estimated 
as the average difference in energy between simple sad- 
dles {K = 1) and threshold minima {K — 0). This 
difference can be extracted from the slope of k{u), as 
AC/(uth) ~ l/[3A:'(uth)] ii]- Note that the slope of 
k{u) for A^ = 70 and A^ = 140 is the same. This is con- 
sistent with the fact that activated processes are local 
in space (as we will discuss later) and therefore barri- 
ers do not depend on the size of the system. We find 
AC/(uth) ~ 2.2 w lOfesTth- This is an important result: 
at the temperature Tth where activation becomes dom- 
inant, potential energy barriers are already very large 
compared to the available thermal energy. Activation is 
therefore highly inefficient at the temperature where for 
the first time it is actually needed. This, we believe, is 
the most striking feature of very fragile liquids and it 
confirms the conjecture of B that the fragility of a sys- 
tem is higher the larger the potential energy barriers at 
Tth- We thus predict that a sharp slowing down of the 
dynamics must occur at Tth- Note that the change in 
the mechanism of diffusion at Tth would not be accompa- 
nied by a dynamic slowing down if barriers at Tth were 
small. In this case there would rather be a fragile-to- 
strong crossover at Tth, as discussed in Q. 

To test our prediction about the slowing down we must 
find the temperature marking the onset of glassiness. For 
soft-spheres, this is generally accepted to be Tc « 0.226 
pd| . This value is affected by the same arbitrariness as 
the experimental Tg, since an arbitrary time scale (set 
by the simulation times) is involved [ pd| . However, as 
stressed in the introduction, the slowing down of fragile 



glasses is so sharp that it makes sense to define aTc, as 
long as one keeps the above proviso in mind. 
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FIG. 2. Van Hove self correlation functions for particle type 
1, at times t = 88 (full line), 177 (dotted), 265 (dashed) and 
353 (dash-dotted). iV = 70. 

Given that we use a non-standard cut-off for the po- 
tential, we perform an independent determination of Tc 
for N — 70. To this end we compute the van Hove self- 
correlation function G" (r, t) from configurations sampled 
in a molecular dynamics (MD) run which uses equilib- 
rium MC configurations as starting points. Of course, 
the MD simulation falls out of equilibrium at higher tem- 
peratures than the MC swap dynamics. G" is defined as 



Gi")(r,i) = — ^(5[r,(t)-r,(0)-r]) 



(2) 



The probability that a particle of type a has moved a dis- 
tance r in a time t is proportional to r'^Gs" {r,t), which 
is plotted in Fig. for several times and temperatures. 

A reliable dynamic diagnostic for Tc is to look at the 
evolution of the first peak of r^Gs (r, t) ||ll[. In the 
liquid phase, the peak moves to the right and rapidly 
becomes Gaussian (top panel of Fig. 0). On the other 
hand, in the glassy phase it takes a huge time to reach 
the hydrodynamic limit, and the simulations show an 
unmoving peak whose area very slowly decreases as a 
secondary peak grows (bottom panel) . The middle panel 
shows an intermediate situation. On this basis, we esti- 
mate Tc w 0.24, not far from the standard Tc. This value 
is consistent with the topological transition temperature 
Tth ~ 0.242 we found above, a result which strongly sup- 
ports our scenario. Let us stress the difference between Tc 



and Tth: The first depends on the time scale of the exper- 
iment and it can be sensibly defined only if the dynamic 
crossover is sharp. The latter marks the point where ac- 
tivation starts ruling the dynamics, and it is uniquely 
defined. In fact, Tth has the same nature as the critical 
temperature of mode coupling theory |17| . 

Our estimate of the barriers can be criticized, since the 
average distance in energy between minima and simple 
saddles neglects the requirement that they must be con- 
nected to each other. To test our estimate we perform for 
iV = 70 a direct sampling of the potential energy barri- 
ers. Starting from a simple saddle we follow the gradient 
in the two opposite directions along the unstable eigen- 
vector, obtaining two connected minima. 
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FIG. 3. Average potential energy barriers as a function of 
the potential energy density of the adjacent minimum. Points 
are an average over 1652 barriers. N = 70. 



In Fig. g|we plot the average barrier size AC/ as a func- 
tion of the energy density u of the adjacent minimum. 
On the same plot we report the value of Wth {N — 70), 
and the estimate of Af7(uth) obtained from the slope of 
k{u): this estimate agrees with the value that can be 
read off from the plot. We conclude that the function 
k{u) provides the threshold energy and the potential en- 
ergy barriers at the threshold. 

A second important piece of information is contained in 
Fig. ]^: the typical barriers grow with decreasing energy. 
The consequence is that below Tth the dynamic slowing 
down is enhanced not only by the decrease of the thermal 
energy available for activation, but also by the increase 
of the typical barriers. Thus, below Th, we expect super- 
Arrhenius behaviour of the relaxation time |lq| . 

A brief comment on barrier crossing is in order here. 
Thermal activation must be introduced within a canoni- 
cal description, since in the micro-canonical ensemble the 
total energy is conserved. In other words, an activated 
transition is performed by a sub-system, with the rest 



acting as a thermal bath. The sub-system must be much 
smaller than the total system, and indeed activated pro- 
cesses involve a finite number of particles. For this rea- 
son, potential energy barriers associated with such pro- 
cesses are finite in the thermodynamic limit, implying 
that simple saddles and minima have the same potential 
energy density for N —^ oo. Yet, the equilibrium poten- 
tial energy density is of order fc^T above minima. This 
fact does not imply that barriers are easy to overcome, 
nor that thermal activation is irrelevant, but simply that 
activated processes involve a finite number of particles. 
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FIG. 4. Distribution of particle displacements. Full line: 
distance between the two minima. Dashed (dotted) line: dis- 
tance between the saddle and the right (left) minimum. Inset: 
same plot in a log-log scale. Dashed-dotted line: distribution 
of largest displacement. A*' = 70. 

As we have shown, activated processes are crucial be- 
low Tth and therefore a real-space description of barrier 
crossing is very important. To this end we computed the 
distribution of the displacement, that is the distance be- 
tween the position of a particle in a minimum and its po- 
sition in the crossing-connected minimum (Fig. m p9[ . 
To interpret this result we need first to fix a reference 
distance: Fig. shows that, even in the glassy phase, 
particles can easily travel a distance r^ ~ 0.5. The pri- 
mary peak in the displacement distribution indicates that 
the large majority of particles moves less than r™ , while a 
small secondary peak can be seen at r « 1 , involving only 
« 2 particles EQ]. From the radial distribution function 
(not shown) we know that r w 1 is the nearest-neighbor 
distance for type 1 particles. These facts thus suggest 
that in this system activated processes involve a small 
number of nearest neighbor particles exchanging posi- 
tions, while many particles move a small amount to make 
way for them. We also measure the largest displacement 
for each pair of minima and find that its distribution 
has no secondary peak at short distances, confirming the 
above interpretation. Finally, the distribution of the dis- 



placements between minimum and intermediate saddle 
shows no secondary peak at large distances, consistently 
with the fact that on the saddle the exchanging particles 
have not completed their transition yet. 

In this Letter we argued that glassy slowing down in 
fragile liquids is caused by the presence of a topological 
transition. Potential energy barriers are much larger than 
the available thermal energy at the transition, and they 
increase with decreasing energy. Activated processes in- 
volve small numbers of particles, each moving a distance 
of the order of the nearest-neighbor separation. 
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